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Abstract. We present the outline of a new model for the 
coherent radio emission of pulsars that succeeds in repro- 
ducing the energetics and brightness temperatures of the 
observed radio emission from the observationally deduced 
distances of 50-100 pulsar radii above the neutron star in 
a narrow region. 

The restrictions imposed by energy conservation, 
plasma dynamics of the coherent radiation process and 
propagation effects are used to apply the action of a 
plasma process like coherent inverse Compton scattering 
(CICS) (see Benford, 1992). In accordance with our find- 
ings (Kunzl et al. 1998a) this process requires Lorentz fac- 
tors of about 10 which are lower than in most other radio 
emission models. This implies that no significant pair pro- 
duction can take place near the surface and we expect 
charge densities close to the Goldreich- Julian value (Gol- 
dreich & Julian (1969)). 

To fulfill the energetic and electrodynamic constraints 
the model requires constant re-acceleration in dissipation 
regions which can be interpreted as a voltage drop similar 
to that in a resistive current circuit built of a battery, 
connecting copper wires and a resistive load. 

Using the emission heights for PSR 0329+54 published 
by Mitra and Rankin (2002) and the spectral data from 
Malofeev et al. (1994) we find that a constant depth of 
the dissipation region of about 2rNS can account for the 
observed luminosities and spectral behavior. 

The extremely high peak fluxes in the substructures 
of single pulses can be explained by beaming effects as 
discussed in Kunzl et al. (1998b). 



1. Introduction 

Pulsars are generally accepted to be rotating neutron stars 
having 10 4 — 10 8 T magnetic fields that induce electric po- 
tential drops along the magnetic field lines which allow 
efficient particle acceleration of thermally emitted elec- 
trons in the region of open field lines (Goldreich and Ju- 
lian (1969); Sturrock (1971); Ruderman and Sutherland 
(1975), Arons (1981)). Although the emitted radio lumi- 
nosity in comparison to the spin-down power of a rotating 



neutron star is tiny (only the 10~ 5 - 10 -6 part of the spin- 
down power is emitted in the radio regime) one of the 
most striking features about pulsars is their very intense 
pulsed radio emission, the physical origin of which is still 
poorly understood. The wealth of observational data (in- 
tensity, polarization, shape of the radio spectrum) clearly 
indicates that the pulsar radio emission cannot be of inco- 
herent origin, but must be due to some coherent radiation 
mechanism (e.g. Lync and Graham-Smith 1998; Kramer, 
Wex and Wielebinski (2000) and references therein). Es- 
pecially at low frequencies (below 1 GHz) a very high de- 
gree of coherence is needed to explain the observed in- 
tensities in terms of brightness temperatures. Therefore a 
pulsar emission model has to explain how a large number 
of phase coupled particles can be stimulated to radiate in 
a coherent manner. 

Currently the most favored models for radio emission 
mechanism are curvature radiation (e.g. Rankin (1992); 
Radhakrishnan and Rankin (1990); Gil (1992)) or maser 
process due to strong turbulence, like the free elec- 
tron maser, (hereafter FEM) (e.g. Melrose (1978), Rowc 
(1992a,b)). 

Since there are severe difficulties with coherent curva- 
ture radiation (Lesch et al. (1998) pointed out that the 
energetic requirements can hardly be matched, Melrose 
(1992) discussed the absence of an appropriate bunching 
mechanism to build up a coherent curvature radiation) 
we propose a process directly coupled to the most signif- 
icant eigenfrequency of the magnctospheric plasma: the 
electron plasma frequency. It has been shown by Melrose 
(1978) that such an amplified linear acceleration emis- 
sion (ALAE) process could produce sufficient growth rates 
i.e. is fast enough to explain the observations. But it needs 
low Lorentz factors of the emitting particles to work. This 
model requires low energy particles with energies of a few 
MeV. 

Recently, support for such a low energy scenario came 
from various criticisms on the so-called standard model 
with an inner gap where a highly relativistic pair plasma 
is created: 

It was argued by Kunzl et al. (1998a) and confirmed by 
Melrose and Gedalin (1999, 2000) that the observed radio 



2 



Kunzl et al.: A current circuit model of pulsars 



emission would not escape from a magnetosphcrc filled 
with a dense relativistic pair plasma. Additional problems 
with the production of high energy particles close to the 
pulsar were discussed by Arendt & Eilek (2000), Lesch 
et al. (2000) and Jessner et al. (2001). They commonly 
concluded that the assumption of strong electric fields at 
the pulsar surface that ultimately result in ultrarelativistic 
particle energies is at least questionable. 

Thus, we introduce an entirely different approach here, 
which has first been promoted by Shibata (1991) for 
an analysis of the global magnctospheric structure of 
rotation-powered neutron stars. In contrast to the stan- 
dard model assuming a shielding of the electric fields by 
massive pair production we propose that in the inner re- 
gion of the magnetosphcrc (i.e. the open field lines inward 
of the neutral surface) the situation is similar to the mul- 
tiple scatterings in a solid state conductor. In other words 
we have a current whose time average is the same at all 
cross sections through the open field line region. There 
is neither acceleration close to the surface (as indicated 
by our recent findings (Jessner et al. 2001) nor deceler- 
ation in the dissipation region (as indicated by the very 
small dissipation rate given by small the radio luminosity 
in comparison to the total spin down power). 

In this conducting wire model any dissipation region 
corresponds to a resistance where electric potentials are 
present. The "deceleration" caused by the losses is bal- 
anced by the "acceleration" due to the electric field. The 
most important consequence of this scenario is that the 
energy loss per particle is no longer limited by the av- 
erage kinetic energy of the electron in the same way as 
the kinetic energy of the electrons in a copper wire cur- 
rent circuit is no indicator for the dissipated energy. The 
one-dimensional approach with a constant voltage applies, 
because the resistivity perpendicular to the magnetic field 
is small inside the pulsar (to allow sufficient supply of 
charged particles for a Goldreich- Julian DC current) and 
large enough in the magnetosphcrc (to avoid "leakage", 
meaning current closure across the field lines in the inner 
magnetosphere) (see Ewart et al. 1975; Itoh 1975). 

Our aim is now to obtain the constraints on the dis- 
sipation field in such a model and see whether we can 
explain the luminosities in the radio frequency range at 
the emission heights given by the observations from a thin 
radiating layer. 

2. The model 

In a current-carrying collisionlcss plasma waves are known 
to introduce resistivity via fluctuating electric fields. These 
fields present the deviations from the ideal plasma state 
with an infinitely high electrical conductivity. For an in- 
coming current I the locations of a significant high plasma 
wave intensity act as a resistance R. The emitted power 
of such a resistive region is given by L — I 2 R. In a time 
averaged picture the current has to be divergence- free. So 



if a DC current flows through a region in which plasma 
waves generate a resistance the particles lose their kinetic 
energy meaning that their velocities decline. Although for 
relativistic particles this effect is small it cannot be en- 
tirely neglected. 

In a dense pair plasma this energy loss of the particles 
presents no difficulties since the velocities of the two com- 
ponents can adjust in an adequate way to fulfill both en- 
ergy and current conservation. However, in a single charge 
plasma deceleration has to be accompanied by a density 
rise. So if the particle beam would gradually lose energy 
the relative density would have to increase all the time. 
This is totally incompatible with the outer boundary con- 
dition of a Goldrcich-Julian magnetosphcrc. 

The latter point is the main difference to the argu- 
mentation by Shibata (1997) who also assumes a current 
outflow determined by the system. In Shibata's model cur- 
vature effects causes space charges to accumulate that can- 
not be removed by acceleration on some field lines as the 
particle velocity is limited by c. We assume that these 
curvature effects do not play a significant role as they will 
be balanced by particles of the opposite charge drifting 
inwards from the neutral line and therefore only deceler- 
ation by losses will result in a space charge. This would 
mean that radiation losses produce a global deviation from 
the GJ case if there was no acceleration of the particles. 
So such a description would immediately turn out to be 
inconsistent (non-idealities induce electric fields which ac- 
celerates or decelerates particles). 

Nevertheless there is a consistent description of a single 
charge plasma distribution if we recall that the rotating 
neutron star induces the vacuum electric field. The elec- 
tric charges on the neutron star surface are displaced in 
such a way that the field that has to be shielded by the GJ 
charge and the field is actually anchored in the neutron 
star. Thus, the vacuum field can be interpreted as a sta- 
tionary distribution of infinitely heavy positive particles 
whose strength is given by the entire pulsar system. 

This means that electrons that suffer losses will imme- 
diately be dragged on by the electric field that builds up in 
such a way that just an electrostatic wave is produced. A 
particle transversing this wave is accelerated by this non- 
vanishing potential drop and looses the gained energy by 
non-thermal radiation. 

In this case the particles have to pile up which means 
that their velocity has to decrease. Where the density goes 
down, the velocity has to increase accordingly to keep the 
current stable (divj = 0). 

Then the resulting electric field is able to accelerate 
particles and varies with half the wavelength of the density 
variation (as follows directly from the Poisson equation). 

Such a system is similar to a battery whose poles are 
the region of the open field lines (acting as the negative 
pole in the parallel rotator) and the edge of the polar cap 
being the anode. If a particle starts close to the pole and 
returns on the last closed field line it encounters a voltage 
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drop of typically 10 13 V. But, of course, this does not 
mean that the particle has to be accelerated by the full 
potential in one step anywhere along its path. 

As will be motivated in sect. 7 it is likely that the 
outer magnetosphere itself provides a constant minimum 
resistance which clearly dominates all plasma induced re- 
sistances so that the latter are only minor corrections. In 
other words, if the main voltage drop occurs close to the 
light cylinder almost all of the energy is dissipated there 
and transformed into high energetic radiation and ultrarel- 
ativistic particles. A schematic view of the model is shown 
in fig. @.. 




Fig. 1. Schematic view of the current circuit model. The 
current is determined by the applied voltage AV and the 
resistances R\ to Ra- Ri is the battery's resistance which 
is assumed to be small. i?2 stands for plasma resistance 
effects, mainly the losses by radio emission. The vastly 
dominating resistance is i?3 which limits the current to a 
relativistic GJ current over the two polar caps. In R4 we 
combine other possible resistances which are also assumed 
to be small (e.g. curvature radiation losses or similar ef- 
fects) 

Of course, there are two basic assumptions in this 
model, which have to be checked carefully. On the one 
hand, the pulsar is described as a voltage source. This de- 
scription requires that the generator's internal resistance 
is negligible compared to the total resistance. It will be 
shown in section 6 that the maximum current allowed 
would be a highly relativistic Goldrcich- Julian current 
over the entire two polar caps. The corresponding resis- 
tance is about 5 Q which is huge in comparison with the 
transverse resistance of about 10 -8 H (Itoh, 1975) Thus 
the first assumption appears to be quite reasonable. 

The second assumption, namely the one-dimensional 
current flow in the inner magnetosphere is also very natu- 
ral. The conductivity can be expressed by a mean collision 
time r using the formula 



The transverse collision time r can be estimated by 
l/w c (the inverse gyro frequency) whereas the collision 
time for the parallel current cannot exceed the inverse 
plasma frequency l/w p (see below). As in the inner magne- 
tosphere uj p <C oj c , the one-dimensional approach is appro- 
priate. Another way of confirming this is to calculate the 
maximum (transverse) drift velocity in the inner magne- 
tosphere. Taking the maximum plasma wave electric field 
(see below) and typical pulsar parameters (P = 0.5 s, 
B = 10 8 T) we find a drift velocity of 

Wdrfft = 0.359 ^x^V 72 (2) 

which is far below c for typical emission heights and 
particle energies used in the model (a; em ~ 50, 7 = 10). 
Particles move almost perfectly along the field lines as long 
as relativistic parallel outflow is possible (on the neutral 
line the situation is a bit different as a pure electron flow 
cannot cross this border easily). 

An important consequence of our description is that 
the total energy a particle radiates along its path is in 
no way limited by its own kinetic energy as the actually 
dissipated energy comes from the inductor, i. e. the bat- 
tery. Here, the analogy to a battery with a bulb attached 
through copper wires can be helpful. In a metal conduc- 
tor the current carrying electrons have velocities of only 
some millimeters per second which means a kinetic en- 
ergy of only 10 _16 eV because of the small path length of 
electrons in a solid state conductor. Nevertheless each of 
them converts a total of 1.5 eV into heat on its way from 
the cathode to the anode. The only restriction which is 
given by the kinetic energy is the potential difference of 
the electrostatic field in one elementary radiation cell. As 
we will point out in sect. 4 this limit directly leads to the 
maximum possible resistivity in the plasma and therefore 
is never exceeded anyway. 

The accelerating fields occurring in the outer magne- 
tosphere do not invalidate the current circuit model, how- 
ever as the circuit is open (high resistance) near the light 
cylinder. 

Another advantage of this description is that our 
model uses particle inertia as the source for an acceler- 
ating field near the pulsar surface to get relativistic par- 
ticles. So a very small acceleration region has to occur 
as the electrons have non-relativistic velocities inside the 
neutron star; therefore an overdense thin charge layer de- 
velops just above the surface where the electrons reach 
their terminal energy. 

Therefore is no need for a "gap" (meaning a starved 
region near the neutron star) but inertia naturally pro- 
vides the fields necessary to reach the mildly relativistic 
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energies, and the terminal energy is determined by the 
global current instead of some large-scale space charge 
field (cf. sect. 6). 

3. Application to radio emission 

In the following we derive a quantitative model for the 
radio emission of pulsars. The particular aim is to explain 
the energetics of the radiation, from observed or observa- 
tionally derived pulsar parameters such as 

— the period P 

— the dipolar surface magnetic field B 

— the brightness temperature Tb of the low frequency 
emission 

— the emission height x em 

— the total luminosity in the radio band L 

With these input parameters we derive the radial 
width of the emission region, which should be small to 
be consistent with observations. A narrow radiation zone 
is suggested by both the highly structured emission in 
space and time as seen by many microstructure experi- 
ments (Boriakoff, 1992). 

3.1. Plasma beams and coherent emission 

Emission processes in relativistic plasmas are always 
connected with the electron plasma frequency, which is 
(Kunzl et al. 1998a; Melrose and Gedalin 1999) 

/ ne 2 
y jeoirie 

As a very promising candidate for the origin of coher- 
ent radio emission of pulsars we consider the collective 
inverse Compton mechanism, which was detected first in 
a number of laboratory experiments in which relativis- 
tic electron beams with a number density rib penetrate 
into a background plasma with density n p and in which 
extremely intense coherent microwave emission has been 
observed (Kato et al. 1983; Levron et al. 1987; Yoshikawa 
et al. 1993; 1994). The detected radiation has been ex- 
plained in terms of collective inverse Compton scattering 
of strong Langmuir turbulence (Benford and Wcatherhall 
1992; Benford 1992). Relativistic electron beams interact 
with self-excited strong density fluctuations in the back- 
ground plasma if the beam-plasma-density ratio rib / n p ex- 
ceeds 0.1 for a beam Lorentz factor of about 3. The beam 
electrons are bunched by the electric field of the density 
fluctuations and radiate coherently in their reference sys- 
tem at jLUpc, transformed into the observer's frame this 
gives an emitted frequency of about j 2 lo pc . The efficiency 
of this energy transfer of beam energy into photon energy 
is up to 30%. 

The general scaling law from laboratory experiments 
for the relation between beam Lorentz factor and density 
ratio is (Benford 1992) 



As shown in Fig. ||for large beam energies the required 
density ratios is quite small, whereas for small beam en- 
ergies the beam density must be comparable to that of 
the background plasma to drive collective inverse Comp- 
ton scattering. For a given beam energy one can calculate 
the required density ratio. 

-Iqqqq . . . . critical Lorentz factor for CICS 



1000 




nb/np 

Fig. 2. Minimum Lorentz factor that allows coherent in- 
verse Compton scattering (CICS) as a function of the ratio 
of beam density rib to background plasma density n p \ (see 
Benford, 1992). For strong beams (rib ~ n p \) even low rel- 
ativistic particles can cause coherent radiation (region I). 
In region II no coherent emission is possible. 

The underlying physics of this mechanism has been 
investigated in all details and the general consensus is, 
that a relativistic electron beam with sufficient energy and 
particle density penetrates into a background plasma the 
collective inverse Compton scattering is the unavoidable 
coherent radiative power output (Benford 1992 and ref- 
erences therein). Based on these findings we transfer this 
process into a neutron star magnctosphcrc and investigate, 
whether this process matches the required properties. 

Collective inverse Compton scattering is a nonlinear 
scattering processes whose radiative output in first order 
does not depend on the origin of the density fluctuations. 
In contrast, direct excitation mechanisms like mascr or 
beam radiation depend sensitively on the form of the beam 
energy distribution and the time evolution of the beam 
(Ursov & Usov 1988; Melrose & Gedalin 1999; Melrose et 
al. 2001). 

For comparison with alternative processes, we give a 
general scaling of the emitted frequency with the nonrcl- 
ativistic plasma frequency uj p °J for the different radiation 
mechanisms 
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(o) 

^pc all 



(3) 



In this notation the value of a for the collective inverse 
Compton scattering is 3, whereas for direct emission it is 
1. In the next chapter we will debate an argument why 
the description of the pulsar system as a resistive current 
circuit makes a = 1 less likely than a higher value. Insert- 
ing the expressions for the Goldreich-Julian density and 
the plasma frequency, we can solve this equation for 7: 
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bJp (^em) 


-2/a ^ 


2ttv 





-l/c 



(4) 



According to the usual estimates, an incoherent pro- 
cess from an electron with a Lorentz factor of 7 cannot 
produce radiation with a brightness temperature exceed- 
ing 10 10 I<7. However, in the case discussed here, this ar- 
gument is not applicable, because there is no absorption. 
The emitted power of a radiation process is not limited 
any more if the radiation produced further inside is not 
absorbed on its way out of the magnetosphere. 

Consider a cylindrical plasma column of radius i?, 
length I, particle density n and single particle power out- 
put Psingic in forward direction. Then the total power adds 
up to 



P gcs = lR z -KnP f 



single- 



(5) 



Thus the brightness temperature can reach arbitrarily 
high values on the condition that I is large enough no 
matter, whether the emission process is coherent or not. 

Consequently we use a different estimate for the bright- 
ness temperature. The average flux is received from an 
area that corresponds to the polar cap moved and ex- 
panded along the field lines up to the emission height: 



Note that there is no additional geometric factor as 
the estimate of the total radio luminosity already averages 
over some typical pulse profile. The values for the bright- 
ness temperature found in (|^) are consistent with observa- 
tionally deduced results (Sutherland, 1979, Kramer, 1995) 
as the same principal method has been used in both cases. 

High coherence is required for a different reason, 
namely that the radiating volume is small enough. In the 
next step we show that our model is able to reproduce 
these fluxes under the assumption of a thin layer emit- 
ting highly coherent radiation, but where the coherence 
factor is kept below the limit imposed by the model. As 
already mentioned, we propose a mechanism of the rela- 
tivistic plasma emission type. 

For the discussed process, particles are scattered by 
strong, nonlinear plasma waves (solitons) excited by the 
two-stream instability. Such a radiation mechanism can 
be described as an inverse Compton scattering (ICS) of 
solitons by relativistic electrons. 

For each particle coherent ICS radiates a power of 



^wavc £ „,2 



(9) 



with the wave electric field E wave if Hu wavc <C m e c 2 , 
which is the case for radio emission^]. N is the coherence 
number, and 7 2 comes from the Lorentz transfomation of 
-E wavc ■ As long as the wave electric field is small, the wave 
can be treated linearly as a small density fluctuation. But 
in an unstable situation the wave grows exponentially. The 
growth saturates quickly because the wave will also accel- 
erate particles as a back reaction. Thus, the wave growth 
has to be calculated using non-linear methods (Benford & 
Weatherall, 1992). As a result, the strongest possible wave 
has an electrostatic field energy comparable to the kinetic 
energy of the plasma. In the extreme, these two energy 
densities are equal 



A = A r 



(6) 



The total radio luminosity over a frequency band of Av 
around the frequency v can therefore be equated to the 
thermal flux of a blackbody with the appropriate bright- 
ness temperature. For broadband emission [y = Av) and 
using the definition of the brightness temperature we find 



L = F v AvA, 



cap J 



k B T B A cap Q 



(7) 



which yields a brightness temperature of 



Lc 2 



1 



;e E 



wave 



which means (j^j can be rewritten as 



rjcoh 
*rad 



NcrTCwy m e c 



(10) 



(11) 



As particles are not assumed to change their longi- 
tudinal momentum in the radio emission region, we can 
interpret the radiated power as an electric dissipation field 
(which is the field necessary to balance radiation losses). 
Its strength is 



E, 



diss 



2irv 3 k B A 
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pcoh 
£rad 

ec 



2Na-T^B eom e c 2 



-7 



(12) 



cm 



1.87- 10 23 K 



cap^em 

L 



10 20 W J V 400 MHz 



■50 



1 For higher wave energies the cross section is no longer de- 
" 3 scribed by ot, so that the Klein-Nishina cross section (ctkn ~ 



or (In "/)/"/) must be used. 
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3.2. Extension of the radiation region and the maximum 
luminosity 

Now we calculate the radial width of the radiation zone. 
This is an important test of the model, as the very small 
RFM taken from the observations requires a fairly narrow 
emission zone. Apart from that, the highly modulated ra- 
diation observed in the radio band would be hard to ex- 
plain for large widths. 

In ( |l2| ) the only the free parameter left is N. However, 
this number cannot directly be fixed from the brightness 
temperature as mentioned before. Nevertheless there are 
upper limits for the coherence factor which serve to derive 
a minimum radial width of the emission region. 

As we consider a plasma wave being the cause of coher- 
ence, the maximum volume of coherently radiating parti- 
cles is a sphere of around one Debye volume. Taking into 
account relativistic effects, Melrose (1992) finds a maxi- 
mum coherence volume of 



3 2 

Vcob. = 1 — 1 — 



(13) 



and therefore 



«GJ — 

.Wo/ 7T 



6.41 • 10 13 7 7/2 



2Q 



P 



1/2 



B, 



-1/2 / Xcm \ 

V 50 J 



B e 5 = 

3/2 



(14) 



Another upper limit is found by equating the parti- 
cle and the field energy density, since in the optimal case 
those quantities are comparable. More detailed calcula- 
tions (Lesch & Schlickeiser, 1987) show that the field en- 
ergy density is less than half of the particle energy density 
(as already used in eq. JlC|)). 

Inserting this result into ( |l2] ) we find a maximum co- 
herence number of 



L 



sjdV 



(16) 



v rad 



where V ra .d denotes the dissipation volume. 

is given by (|l2|) and (|l5|), j(x) = n G j(x)ec, and 
by assuming that the non-ideality extends over the entire 
cross-section of open field lines; using a dipolar geometry 
and neglecting angular corrections to the GJ- density all 
quantities in ( |Tg| ) only depend on the radial coordinate. 
The volume element dV^ can be transformed to 



dV = r NS A cap x 3 dx 



(17) 



Here we have used (g) , as the cross-section is a simple 
function of x. Therefore ( [l6|) reduces to a one-dimensional 
integral. 

Inserting all the previous expressions the integral reads 



27r 2 J/ 2 m e 
Boile 



l/(2c 



Xam+'o 



2ttB m c c 2 2QB e 



-ec 



x3 /(2a) x - 3 /2 da . 



(18) 



where lo denotes the radial extension of the radiation 
region in units of t^ns, whereas x CTn stands for the height 
of its inner edge in pulsar radii. 

Evaluating the integral and inserting numbers for a 
frequency of v = 400 MHz we find 



10 22 W | — 



I/q-5/2 



B 3 8 /2 ~ 1/a F Q (19) 



where 




N 



Bs 1/2 (^-) 3/2 • (15) 

Comparing the maximum coherence numbers from 
( |l4| ) and ( [To] ) we find that the latter is smaller whenever 
7 > 1.97. Therefore it is a good approximation to use ( |l5| ) 
for determining the maximum coherence number. 

For estimates of the radial extension of the radio emis- 
sion region we take the total radio power and calculate 
the radiating volume necessary to reproduce the observed 
luminosities. 

The total energy loss in a resistive current reads 



F a = 0.0286 1/a < 4 



r\l/4 
10 



1/4 



ln^ 



i+lo 



for a 
for a 



for a = 3 

(20) 



Solving this equation for Iq in the case a = 3 renders 

= Xem (21) 



exp 



[ - 



13/6 



B, 



7/6 / L 



Vio 20 w 



- 1 



which, for typical parameters (L = 10 20 W, P = 1 s, 
B 8 = 1, x om = 50) yields l = 0.295 which justifies the 
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radial extension vs. luminosity 
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Fig. 3. Radial extension of the radio dissipation region 
(in units of pulsar radii) versus luminosity for a standard 
pulsar (P = 1 s, Bg = 1). The emission height x cm has 
been set to 50 pulsar radii. 

assumption of a narrow emission region. The dependence 
of and L is shown in fig. ||. 

From the result above one can obtain sufficient lumi- 
nosity even for slow pulsars. However, for small O the ra- 
dial extension of the emission region increases. Especially, 
the corresponding electric field (|l^) can even exceed the 
vacuum field [E vac = QB rwsx~ 4 ) locally. In a shielded 
magnetosphere the kinetic energy density and by equipar- 
tition the dissipated electric field, may exceed the energy 
density of the vacuum field. 

4. Application to PSR 0329+54 

The pulsar PSR 0329 is an extensively researched object 
and one of the few where the spatial orientation of mag- 
netic and spin axis is well determined (Rankin, 1983). As- 
suming that the profile peaks occur on the same field lines 
for all frequencies enabled Mitra and Rankin (2002) to ob- 
tain revised effective opening angles and emission heights 
for PSR 0329+54 amongst others. A good fit to the sur- 
face opening angle (including relativistic corrections) for 
the field lines that provide the two outer peaks of the pro- 
file is (f>o — 0.66°. For any angle <j> between those peaks we 
can find an emission height of a profile at a specific fre- 



quency using the parametric equation of dipole field lines: 

ls sin(</>) 



r{4>) 



/i sin(<£ ) 



(22) 



Mitra and Rankin (2002) found that the emission heights 
range from about 500 km for frequencies above 1 GHz up 
to 1100 km for the lowest frequency of 100 MHz. Together 
with the known radio fluxes (Kramer et al. 1997) we were 
able to match the described model to the data and deter- 
mine the dissipation length scales for frequencies between 
100 MHz and 30 GHz. To enhance the realism of the fit, 
the relativistic (Muslimov & Tsygan, 1992) corrections to 
the Goldreich-Julian densities and effective areas were also 
included. In this case we calculated the local Goldreich- 
Julian density through 



n(r) = 2e B 



(23) 



Here 



T 

-3(f) 



z r 



r 



r 



and fi = /(^ns) are the corrections due to general rel- 
ativity (Muslimov & Tsygan). In this context r g = 2r ™"° 
denotes the Schwarzschild radius of the neutron star. 

Assuming that the observed emission is indeed caused 
by coherent inverse Compton scattering, we can calculate 
the Lorentz factor of the emitting electron current as a 
function of distance from the pulsar. 



7( w obs 



n(?"em(Wobs))e 2 



to 



obs 



(24) 



It varies from 13 (100 MHz) to 267 (30 GHz). Using 
the relativistically corrected values for distance, density 
and Lorentz factor we obtain an improved estimate of the 
dissipation scale A which is given by 



A 



P co h(r, 7) ' n i r ) ■ ^r 2 (1 - cos(0)) 



(25) 



Here L and 4> are the observed luminosity and opening 
angles of the pulsar beam at a given frequency and P co h 
means the emitted power per coherently radiating particle. 
We estimated the luminosity as an upper limit from the 
radio flux by assuming that the smaller cap area given by 
the peak flux opening angle is filled with uniform emission 
can be a good representation of the conal beam structure. 

Fig. 4 shows the development of the the dissipation 
scale with emission height. After a steep rise at 500 km 
{v > 1000 MHz), A stays around 20 km up to a height 1100 
km (100 MHz). Thus we have evidence for a quick onset 
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where P ra( j {v) is the total power emitted in frequencies 
below v , whereas d denotes the transverse extension of 
the radiation region. For a rough quantitative estimate, 
we can set P r ad( I/ crit) = L, and replace the derivative by 
a quotient. Thus we obtain the simpler equation 



Sobst^crit) 



L 



(27) 



Fig. 4. 

height. 



Dissipation length A as a function of emission 



This flux is drastically enhanced by two effects, pro- 
ducing a strongly anisotropic emission pattern. The first 
is the well-known relativistic lighthouse effect (Rybicki 
& Lightman, f979). A relativistically moving source of 
isotropic radiation emits most of the power in a narrow 
cone with an angular opening of about I/7 for 7 3> 1, as 
can easily be verified by performing the Lorentz transfor- 
mation. The exact result is 



of a strong and stable coherence at a height of 600 km, 
lasting for another 600 km. This is also borne out by the 
fact that the brightness temperatures of the source vary 
only a little in this regime where most of the pulsar ra- 
dio luminosity is created. A typical fluctuation time scale 
of 7 • A/c = lms is associated with the dissipation scale 
below 1 GHz. Above that frequency it is expected to de- 
crease monotonically to 20/iS at 20 GHz. These timescales 
are in good agreement with observational evidence (Langc 
et al. 1998) where typical microstructure timescales were 
found to range within 600 — 1200/is for that pulsar at 
4.85 GHz. There is also some observational evidence for a 
decrease in microstructure timescales with increasing fre- 
quency for the pulsars PSR B1133+16 and PSR B2016+28 
(Lange et al. 1998). 

5. Micropulse flux enhancement by beaming 
effects 

In this section we discuss anisotropy effects for pulse sub- 
structures which may enhance the observed brightness 
temperatures by up to several orders of magnitude. It is 
very unlikely that the flux enhancement only arises from 
a growth of the dissipation region. For the nanosecond 
giant pulses of the Crab pulsar this is even impossible 
(if the width of the radio emission region jumped from 
only 0.001 — which is extremely narrow — to around 1 
pulsar radius, the peak would occur on a timescale of 
T m j C W hms/c ~ 3 • 10 -5 s). But as such small time scales 
also mean highly localized emission, the anisotropy of the 
elementary radiation process has to be considered. For an 
isotropic process there is a simple connection between the 
radiated power and the observed flux: 



(26) 



arctan 



sin# 



7 (cos6> + /3) 



(28) 



with 9 being the angle between the direction of the 
flux and source velocity in the co-moving frame and 9' 
denoting this angle in the observer's system. For 9 = 90° 
and 7> 1 the approximation 9' « I/7 is a good one. 
Therefore we find a flux enhancement of 
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47T 

tt/7 2 



= AY 



(29) 



by the lighthouse effect. 

Another anisotropy factor is the coherent radiation 
process itself, as has been shown by Kunzl et al. (1998b). 
To understand this effect, one has to recall that coherent 
emission requires a phase coupling of the emitting par- 
ticles in one direction. For relativistic plasma emission, 
the preferred direction obviously is the direction of the 
streaming velocity (forward direction). 

As long as the spatial dimension of the coherently radi- 
ating volume is small compared to the emitted wavelength, 
the interference is still constructive even under large angles 
to the forward direction. However, as soon as the exten- 
sion of the coherence region becomes comparable to the 
wavelength, the coherent emitter can be seen as a phased 
antenna field producing intrinsically beamed emission. 

The quantitative analysis of this effect has been per- 
formed by Kunzl et al. (1998b). The first minimum ap- 
pears under the angle 



1 



2tt 



1 



(30) 



Here T := 2ird/ A with d being the transverse extension 
of the emission region. The angle of the minimum is much 
smaller than unity for T ^> 1. 
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Similar to the lighthouse effect, we can use this an- 
gular extension to estimate the enhancement factor by 
coherence beaming. Thus the flux grows by a factor of 



4-7T 



1+^-1 



(31) 



A more realistic treatment of the geometry uses a 
cylindrically symmetric soliton with a density profile pro- 
portional to cosh _2 (i?o/i?) where Rq ~ Ad is the typi- 
cal length scale. Numerical calculations for this structure 
show a much stronger anisotropy (Kunzl et al. 1998b) 

Combining the two anisotropy effects we find the total 
beaming factor 




(32) 



which, even with quite moderate parameters (like 7 = 
10,x cm = 50, N = 10 13 for a typical pulsar), causes an 
enhancement of some 10 3 . Therefore it can easily explain 
strong substructures, and micropulse brightness tempera- 
tures that are some 10 3 times above the mean value (Bo- 
riakoff, 1992). 

As these beaming effects apply only to the elementary 
emission process, but do not enhance the mean flux sig- 
nificantly, it is an important check on whether the model 
can reproduce the observed brightness temperatures of up 
to 10 31 K (Hankins, 1996) observed in giant pulses of the 
Crab pulsar. 

Therefore we take the maximum possible number of co- 
herently radiating particles, and the corresponding beam- 
ing effects. Let A be an arbitrary cross section. Then 
starting with eq. (||) one can express the integrated flux 
I := I^Av = i ra dio,A/^l by the brightness temperature: 



L 



radio, A 

~A 



2ttu 3 k B T B 



(33) 



^crit 



Inserting the critical frequency of the Crab pulsar v = 
160 MHz and T B = 10 31 K, the integrated flux is 



/ = 3.95-10 16 ^ 
nr 



(34) 



As this value means the beamed flux, the actual power 
per area is only 



I ,K W , 

/real :=- = 2.47-10 15 — 7 ~ 2 
£ m^ 




(35) 



Now we can compute the minimum thickness of a layer 
producing the power per cross-section. To obtain the lower 
limit we take the strongest dissipation field possible (see 



eq. (|12j) ) and assume a coherence cell with a lateral exten- 
sion of c/i/p, so that T = 2^7. 

The total dissipated power from a volume with the 
cross-section A and the radial extension d is 



-Pdiss — ^real^ — -Ediss dj A. 



(36) 



For a relativistic Goldreich-Julian current j = jc,.i = 
no]{x)ec eq. ( [36] ) can be solved for d. Inserting the Crab 
pulsar parameters (P = 33.4 ms, B = 3.8 ■ 10 8 T, v = 160 
MHz, x em = 80) and additionally using eq. (Q) we obtain 



d = 3.85i 



( r B 



V10 31 K, 

which corresponds to a time scale of 



(37) 



1.28- 10 _s s 



10 31 K 



12.8 ns 



Tn 



10 31 K 



(38) 



Thus we can expect to see giant pulses on nanosecond 
timescalcs, with brightness temperatures above 10 30 K, 
from far inside the light cylinder in the Crab pulsar, al- 
though the relative emission height will be considerably 
larger than for average pulsars. 

6. Densities in the radiation region 

Here we briefly discuss if the model described above could 
also work for particle densities significantly above or be- 
low the GJ value. We will show that even if there were a 
mechanism producing a current density and particle ener- 
gies matching the requirements to produce the low radio 
frequencies, such a scenario would be very unlikely. 

Assuming that the particle density was much larger 
than the value predicted by Goldrcich & Julian (1969), 
the Lorentz factor corresponding to a certain frequency 
would drop according to (^) . Especially for young pulsars 
this would quickly require 7 < 1 which is, of course, not 
possible (even (3j > 1 should be fulfilled as the particles 
are definitely relativistic). This important fact has already 
been pointed out by Kunzl et al. (1998a). 

A relativistic current significantly below the GJ cur- 
rent (meaning that the particle density is smaller than uqj 
faces even more difficulties. On the one hand, accelerating 
fields are incompletely shielded (the equilibrium Lorentz 
factor where acceleration and radiation losses — now in- 
cluding coherent radio emission — balance is extremely 
high, so coherence is highly doubtful). Furthermore the 
one-dimensional approach would no longer account as a 
permanent strong deviation from ideality allows efficient 
E x B drift. But also the existence of strong 7 -pulsars 
provides an argument against this solution: 

As we observe pulsed 7-emission of pulsars whose high 
energetic luminosity reaches up to 10% of the spindown 



10 



Kunzl et al.: A current circuit model of pulsars 



luminosity, some mechanism must efficiently convert parti- 
cle energy into (hard) radiation (Kanbach, 2001). But the 
potential difference a particle can suffer is limited by the 
polar cap potential. Since the power dissipated by a rela- 
tivistic GJ current flowing over the entire polar cap area 
with the polar cap voltage is comparable to the spindown 
power of the neutron star, a large fraction (well above 
10%) of the generator's power must be pumped into par- 
ticles. This, of course, is only possible for currents of about 
the GJ current and thus for particle densities comparable 
to the Goldreich- Julian value. 

7. Discussion 

We have presented the outline of a model which may 
explain the creation of radio emission within the en- 
ergy constraints of the observationally deduced emission 
heights of 50-100 pulsar radii, (e.g. Blaskiewicz et al. 1991; 
Taylor et al. 1993; Kijak and Gil 1997,1998; Kramer et 
al. 1997). Several other pulsar models fail to explain this 
feature, e.g. the model of Lyutikov et al. (1999) which uses 
Chcrenkov drift resonance, but predicts the radiation zone 
to be in some 1000 pulsar radii. 

One of the main differences to the standard model is 
that we do not assume an inner gap which accelerates 
particles to ultrarelativistic energies and is saturated by 
a pair avalanche. Instead we suggest a conducting wire 
model with a quasi- stationary, low rclativistic current 
that suffers fluctuations on short time scales. Such a model 
is supported by calculations on particle densities available 
from the polar cap by thermal and field emission (Jessner 
et al. 2001) which show that inner gap models have severe 
difficulties in producing particles energetic enough to start 
pair cascades for parallel or nearly parallel rotators. 

We showed that the brightness temperatures for the 
mean pulse can be reproduced by our model and also mi- 
cropulse fluxes can be explained by appropriate beaming 
assumptions. 

The radiation process in the framework of such a model 
could be some relativistic plasma emission like the ALAE 
process described by Melrose (1978). This mechanism re- 
quires mildly relativistic particles to be efficient. As shown 
also in Kunzl et al. (1998a) any plasma process responsi- 
ble for the low frequency emission in the Crab pulsar must 
work with mildly relativistic particles. Von Hoensbroech 
et al. (1998) showed that in a Goldreich- Julian magneto- 
sphere with low relativistic particles propagation effects 
can also qualitatively explain the observed polarization 
properties. 

Finally we would like to note that our model does not 
rule out pair production in general. As can be seen from 
several observations of the Crab pulsar a dense particle 
wind is needed which transports energy to the nebula. The 
calculations of Hirotani and Shibata (1999) show that a 
significant fraction of the vacuum potential drop can occur 
in an outer gap if the X-ray photon density is not too 



high. For a dctaillcd geometric examination of the outer 
gap structure see Cheng et al. (2000). 

It has been shown that for reasonable parameters an 
outer gap model could easily explain the observed energy 
and particle fluxes in the pulsar wind. Recently Crusius 
et al. (2000) have shown that an outer gap model can 
explain the observed spectrum and energy output of the 
Crab pulsar from the infrared to X-ray frequencies which 
provides another hint that the main dissipation region is 
in the outer magnctosphcrc. 

We gratefully acknowledge constructive discussions 
with J. Arons, D.B. Melrose M. Kramer, D. Mitra and A. 
von Hoensbroech. In addition D. Mitra and M. Kramer 
provided us with access to their radio flux- and emission 
height measurements for PSR 0329+54 for which we are 
particularly grateful. 
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